Intrinsic stickiness in open integrable billiards: tiny border effects 
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Rounding border effects at the escape point of open integrable billiards are analyzed via the 
escape times statistics and emission angles. The model is the rectangular billiard and the shape 
of the escape point is assumed to have a semicircular form. Stickiness and self-similar structures 
for the escape times and emission angles are generated inside "backgammon" like stripes of initial 
conditions. These stripes are born at the boundary between two different emission angles but same 
escape times. As the rounding effects increase, backgammon stripes start to overlap and the escape 
times statistics obeys the power law decay and anomalous diffusion is expected. Tiny rounded 
borders (around 0.1% from the whole billiard size) are shown to be sufficient to generate the sticky 
motion, while borders larger than 10% are enough to produce escape times with chaotic decay. 

PACS numbers: 05.45.-a,05.45.Ac 
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I. INTRODUCTION 

Experiments usually measure a signal (light, atoms, 
particles current etc) which comes out from the system 
of interest. If the open set (the escape point or hole) of 
the physical device, where the signal comes out, has a 
relevant size and shape, the signal may include informa- 
tions from inside the physical device and from the open 
set itself. The shape and size of the open set depends on 
experimental interests but also on how the device is build. 
For the dynamics in mesoscopic systems and nanostruc- 
tures, for example, the shape and the size of the small 
open set can present irregularities or defects which may 
induce undesirable changes in the outcoming signal. Such 
sometimes intrinsic irregularities affect the dynamics of 
particles which collide with them. 

From the theoretical point of view it is very difficult to 
describe, in general, the dynamics of colliding particles 
with irregular boundaries. Therefore in recent years more 
and more attention has been given to the description of 
particles confined inside boundaries (or billiards) which 
present some specific edges, softness etc. To mention 
some examples we have the edge roughness in quantum 
dots [i| , unusual boundary conditions in two-dimensional 
billiards Q, Q , effects of soft walls [4] and edge collisions 
M of interacting particles in a ID billiard, rounding edge 
[Sl and edge corrections in a resonator, deformation 
of dielectric cavities [8| , edge diffractions and the corre- 
sponding semiclassical quantization [ol, [IQI , among oth- 
ers. 

Different from the above works, which focus on the 
boundaries of systems, here we analyze the effect of ir- 
regularities from the open set itself. What is the effect 
of (rounded) open sets on the outcoming signal ? To 
mention some examples, rounded open sets are common 
in experiments with semiconductors devices [ll|, [l2| and 
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quantum cavities [l3|, where the open sets have a shape 
very similar to those shown in Fig. [TJ which is the model 
used here (Type I and II). In the above mentioned ex- 
periments the ratio between the open set (not its width, 
but the radius of the rounded border) and the whole de- 
vice lies around ~ 0.001 for Type II borders (the smallest 
case) and around ~ 0.3 for Type I borders (the largest 
case). Since experiments measure a signal which comes 
out from the system of interest, how is it affected by 
such rounded borders ? is it possible that tiny rounded 
borders transform an integrable dynamics into a chaotic 
one ? what is their infiuence on the dynamics in open 
nanodevices, conduction fiuctuations in semiconductors 
[ll|,[l^, particles transport in nanostructures, cold atoms 
in open optical billiards [15| etc? Recent works |16l-[l9| 
in this direction analyzed the effect of the width of the 
open set on the escape rates of particles in open billiards. 

Using the example of a billiard model we show here 
that the shape or structure of the open set induces sticki- 
ness [20|, self-similar and a chaotic output signal, even if 
the dynamics inside the billiard is originally regular. We 
consider that the device borders (or extremities), which 
delimit the holes in realistic open systems, are not sin- 
gle points, but have their own shape. In this work re- 
sults are shown for the escape times (ETs) statistics and 
the emission angles Of in the open rectangular billiard 
shown in Fig. [H While self-similar structures are clearly 
visible in the ETs and escape angles, the sticky motion 
is observed by the power law decay of the ETs statis- 
tic. Although it is possible to detect and quantify sticky 
motion via the distribution of finite time Lyapunov ex- 
ponents [3, in, l2ll - [23| , for the purpose of the present work 
it is more adequate to use the ETs statistics. The fractal 
behavior of the ETs dynamics was also shown recently in 
Bose-Einstein condensates [2^, tr app ed ultracold atoms 
|25j and in a vase-shaped cavity [26|, in generic chaotic 
cavities [23, just to mention some examples. 

The paper is organized as follows. Section ITTl presents 
the model, defines the ETs used to detect the sticky be- 
havior and shows related numerical results. Sections IIIII 



and HVl show the rich dynamics generated by the rounded 
border: long-hved states and self-similar structures for 
the ETs and emission angles. In Section |Vl we present 
our final remarks. 



II. 



ROUNDING BORDERS GENERATING 
STICKINESS 



ETs 



statistic is defined [28| by Q{r) 



The 
limAT^oo ^^ where A^ is the total number of trajecto- 
ries which escape the billiard and A^^- is the number of 
trajectories which escape the billiard after the time r. 
For systems with stickiness the ETs statistic decays as a 
power law [28| Q{r) ex r~^^^^, where jesc > 1 is the scal- 
ing exponent. For hyperbolic chaotic systems and long 
times the ETs statistic decays exponentially. It is known 
that the diffusion exponent /i from the mean square dis- 
placement of the position 
M = 3 -7esc [2i,|30J. 



^^) ^ t^/is related to jesc via 



initial conditions, the dynamics should be equivalent to 
those obtained in the Sinai [3ll and Bunimovich stadium 
|32| billiards. These billiard have marginally unstable 
periodic orbits (MUPOs), which are orbits bouncing per- 
pendicularly between the parallel walls, and are known 
|19l , [33| to generate the sticky exponent jesc = 2 for the 
ETs. In the rectangular billiard the MUPOs are obtained 
from the initial conditions = 0,7r/2. Therefore, in our 
case we also expect to see jesc ^ 2 as R/L increases. 

We start discussing numerically the quantity Q{r) for 
the Type I billiard and for different values of the ra- 
tio R/L, where L is kept fixed. Results are shown 
in Fig. [1 for R/L = 0, 1/10000, 1/1000, 1/100, 1/10 
and 1. Note that these results are automati- 

cally applied to the Type II billiard with r/L ~ 
0,1/200000,1/20000,1/2000,1/200 and 1/20, respec- 
tively. The only difference is that for the type II billiard 
the escape time r' is related to r by r + tp, where tp 
is the time the trajectory needs to travel the two par- 
allel boundaries from the escape hole (see Type II bor- 
der in Fig. [T]). First observation is that for R/L = 0, 




FIG. 1. (Color online) The rectangular billiard, with dimen- 
sion L X D, showing the open sets (borders) of Type I (with 
radius R) and Type II (radius r ^ R/20) considerer here. The 
escape point lies exactly in the middle of the billiard and has 
constant aperture a — Ix 10~^. Initial angle Oo and, schemat- 
ically, the shortest escape trajectories (^i, A2, A4 and Aq) are 
shown. In all simulations we use L = 4 and D = 10. 

In the simulations, particles start at times t = from 
the escape point with an initial angle ^0 towards the in- 
ner part of the billiard with velocity \v\ = 1. The particle 
suffers elastic collisions at the billiard boundaries and at 
the rounded border of the open set. For each initial con- 
dition we wait until the particle leaves the billiard and 
record 6f and the ETs. We use 10^ initial conditions 
distributed uniformly in the interval 0.10 < ^0 ^ 1.54. 
The dynamics for 0.0 < ^0 ^ — 7r/2 is symmetric. With- 
out rounding effects {R = r = 0) the rectangular billiard 
is integrable and has zero Lyapunov exponents. As the 
ratio R/L increases, the rounded border acts like a dis- 
persing boundary generating chaotic motion without re- 
gular islands. For larger values of R/L and some specific 
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FIG. 2. (Color online) Behavior of Q{t) for different values 
of the ratio R/L. 

Q(t) did not has a power law (neither an exponential 
decay). This is the integrable case and the maximal ET 
found is r ~ 5.7 x 10^. Such maximal time is obtained 
because we did not start initial conditions close to the 
MUPOs {0 = 0,7r/2). Separately (not shown) we used 
only initial conditions very close to the MUPOs. We 
found that Q{r) did not have a finite maximal time, but 
a power law tail for very long times with 7esc ^ 1- This 
shows that the MUPOs without the dispersing compo- 
nent (and no chaotic motion), does not generates the ex- 
ponent 7esc "^ 2. This agrees with results obtained [19| 
for some specific initial conditions in an open billiard. 

For very small rounding effects, R/L = 1/10000, the 
qualitative behavior of Q{r) starts to change (when com- 
pared to the case R/L = 0) for times r > 5.7 x 10^, 
i. e. for those trajectories which stay longer inside the 



billiard. This means that the very few trajectories which 
collide with the rounded border tend to stay longer inside 
the billiard and also change the qualitative behavior of 
Q{t). The fitted escape exponent is jesc ~ 0.6. By in- 
creasing the border to R/L = 1/1000, we observe a power 
law decay in Fig. [2] for times r > 5.5 x 10^. The escape 
exponent is jesc ^ 1.3. For R/L = 1/100 and r > 3 x 10^ 
we obtain jesc ~ 1-8 and for R/L = 1/10 and r > 1 x 10^ 
we obtain jesc ~ 2.1. From these points is it possible to 
get 7esc ^ 2.5 + 0.17 ln(i?/L) and consequently we obtain 
the anomalous exponent /i ~ 0.5 + 0.171n(i?/L) for the 
sticky region. 

The ETs from the long living trajectories present 
significant characteristics of sticky motion for R/L = 
1/1000, 1/100, 1/10. In other words, sticky motion and 
long living trajectories start to occur for very small 
rounding borders: around 0.1% from the whole billiard 
size is sufficient to generate the sticky motion and change 
the output signal. Visually such borders are almost neg- 
ligible. Take for example the border in Fig. (TJ it has a 
ratio R/L - 1/143. 



III. ROUNDING BORDERS GENERATING 
THE RICH DYNAMICS 

The physics involved in the dynamics becomes more 
evident when the ETs (r) and escape angles Of are plot- 
ted as function of the initial incoming angle 6^0 and for 
different ratios R/L. These plots are shown in Fig. [3] 
and were generated by using 500 x 500 points in the 
intervals 0.01 < 6>o < 1.0 and 0.00355 < R/L < 1.0 
[-8.0 < \og{R/L) < 0.0]. We start discussing Fig.[3i;a), 
where each color represents a given value of the ETs writ- 
ten as log(r) (See the colobar on the right). Horizontal 
stripes with different colors are evident for a significant 
range of R/L values. Each stripe is defined by a bunch 
of initial conditions which leave to the same ET and con- 
sequently have the same color. For example, for some 
specific initial angles (^o ^ 0.39, 0.56, 0.69, 0.89) we ob- 
serve dark blue stripes which correspond to very short 
ETs. For R/L = these angles can be obtained analyti- 
cally for periodic orbits (for the close billiard) with period 

2n. They are 6>q = arctan [2^] ? where n = 1, 2, 3, . . ., 
which is the middle value of the main dark blue stripes 
which are born at R/L = 0. The corresponding ETs 
are tn = nt, where t is ETs from the case n = 1, ex- 
plained below. The shortest ET {to ~ 8, 0) occurs for 
^0 ^ 0.0, where the particle collides once against the wall 
in front of the escape hole and than leaves the billiard 
[see trajectory Ai in Fig. [l]. The width of the stripes 
are always related to the aperture a from the hole. The 
next shortest ETs {ti = t ^ 12.8) occurs for n = 1 
at O'^^^ = arctan 5.0/4.0 ~ 0.89 where the particle col- 
lides directly with the edge of the billiard. See trajectory 
A2 in Fig. [II The third shortest ETs (- 2t) occurs for 
^(2) = arctan 5.0/8.0 ~ 0.56 and is shown by trajectory 
A4 in Fig. [H We classify these trajectories as A2n and, as 



we will see later, are the starting point from the Arnold 
tongues. In the limit ^0 -^ 7r/2 the ETs tend to increase 
since the trajectories are parallel to the escape point. The 
ETs statistics for these trajectories will not be considered 
here, since the large ET is artificially created due to the 
location of the escape point in the horizontal axis. As n 
increases the ETs from the trajectories A2n increase, the 
corresponding stripes assume other colors (light blue -^ 
green -^ yellow) [see Fig.[3l^a)] and their widths decrease. 
This is the main behavior of the ETs and stripes close to 
R/L - 0.0. 




FIG. 3. (Color online) (a) Log of the escape time and (b) 
emission angle Of as a function of log {R/L) and ^o- 

Increasing R/L we observe in Fig. [3l^a) that the stripes 
for shorter ETs (dark blue) tend to survive longer the 
bounding border effect, while the stripes for larger ETs 
(light blue, green and yellow) tend to be destroyed, or 
even to mix themselves. Notation: we call the dark 
blue stripes as the main stripes and the other stripes 
as the secondary stripes. Therefore, secondary stripes 
are always related to intermediate and larger values of 
ETs. Arnold tongues [34] are visible in Fig. [3fa) for 
log {R/L) > —4.0 and are born around the main stripes 
with frequencies 1/1, 1/2, 1/4, . . .. Their relation with the 
trajectories from Fig. [T] is: Ai -^ 1/1, A2 -^ 1/2, A4 -^ 
1/4, .... Outside the main stripes the dynamics becomes 
very rich and complex. Before explaining how this oc- 
curs, we would like to bring to attention the rich dynam- 
ics generated by the rounding effects. 

Figure ^h) shows Of as a function of log {R/L) and 
Oq. Each color is now related to one emission angle Of 
(see the colorbar). These emission angles vary between 
Of ~ 1.4 (almost horizontally to the left) and Of ~ 4.5 
(almost horizontally to the right). As in Fig. [3l^a), also 



here stripes with different colors are evident for a signif- 
icant range of R/L values. Each stripe is defined by a 
bunch of initial conditions which leave to the same Of. 
In most cases these stripes occur for the same values of 
^0 from Fig. [3l^a). However, two stripes with the same 
color (same ETs) in Fig. [3fa) have not necessarily the 
same color (same Of) in Fig. [3jb). In other words, dif- 
ferent escape angles can have the same ETs. As R/L 
increases more and more, some stripes survive while the 
other ones are destroyed or mixed, as in Fig|3l^a). The 
emission angles show a very rich dynamics due to the 
increasing rounded borders, alternating between all pos- 
sible colors. This will be discussed below in more details. 



IV. ROUNDING BORDERS GENERATING 
BACKGAMMON STRIPES 

In Fig.m^a) and (b) we show, respectively, a magnifica- 
tion for log (r) and Of close to the first region. The mag- 
nification is taken around the trajectory A2 (see Fig. [T]) 
which has a short ETs (main stripe). We observe in 
Fig. [H^a) that above the main stripe the ETs dynamics 
changes significantly when R/L increases: larger (light 
blue and yellow) and shorter (dark blue) ETs appear in- 
side a secondary stripe. This secondary stripe increases 
linearly its width with R/L. For the emission angle [see 
Fig. m^b)] we see a very rich dynamics emerging inside 
such a secondary stripe, alternating between all possi- 
ble colors (all emission angles). In addition, below the 
main stripe a sequence of secondary stripes appear in 
the light blue background, as can observed in the magni- 
fication shown in Fig. [5fa). The width of each stripe in 
this sequence increases with R/L^ remembering stripes 
from a backgammon board. Notation: secondary stripes 
with increasing width will be referred as "backgammon 
stripes" . Stickiness and the power law behavior observed 
in Fig. [2] are born for initial angles which start inside 
the backgammon stripes. These are the initial condi- 
tions which collide, at least once, with the rounded bor- 
der. The emergence of the power law behavior becomes 
more evident if we compare Fig. [5fa) with the emission 
angle behavior shown in Fig. [5fb). The light blue back- 
ground observed in Fig. [5](a), which corresponds to one 
ETs, has two colors (blue and orange) in Fig.[5fb), which 
correspond to two emission angles (~ 2.1 and ~ 3.5). In- 
teresting is that the sequence of backgammon stripes in 
Fig.ISKa), and the corresponding multicolor backgammon 
stripes from Fig. [Sjb), are born exactly at the boundary 
between the blue and orange escape angles at R/L ^ 0. 
Inside the backgammon stripes the range of allowed ETs 
increases very much. This can be observed by noting the 
increasing number of yellow and red points. The dyna- 
mics involved in the emission angles inside the backgam- 
mon stripes is also impressive, showing that tiny changes 
or errors in the initial angle may drastically change the 
emission angle. 

The key observation here is that the dynamics inside 
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FIG. 4. (Color online) Magnifications from the boxes shown 
in Fig. [3] (a) and (b), respectively. 



the backgammon stripes is the consequence of trajecto- 
ries which collide with the inner part of the semicircle 
from the escape point, generating the power law behavior 
for Q{t) (see Fig.[2j). Another observation is that the lo- 
cation of the backgammon stripes itself is not self-similar 
but inside the backgammon stripes the self-similar struc- 
ture is evident. Many simulations (not shown) were per- 
formed to check this statement. As R/L increases more 
and more, the self-similar structures increase very fast, 
always in form of backgammon stripes which emerge at 
different initial angles at R/L ~ 0.0 (this was checked for 
many other initial angles). A simple geometrical calcula- 
tion determines the angle interval [— ^*,^*] for which at 
least one collision with the rounded border occurs, where 



arctan 



smof 



(1 -cos a) ^a/2R 



and a is the collision angle of the trajectory with the 
rounded border, defined relative to the center of the semi- 
circle. The interval [— 6>*,6>*] defines the border lines of 
the backgammon stripes. The larger ETs which appear 
inside the backgammon stripes start to dissapear when 
the chaotic region is reached close to R/L ~ 1.0. 



V. CONCLUSIONS 

Since experiments usually measure a signal coming out 
from the system of interest, the description of open phys- 
ical devices is of utmost relevance. For experiments at 
the frontier of technological limitations, the open set (the 



lO.O 




FIG. 5. (Color online) Magnifications from the boxes shown 
in Fig. m respectively. 



escape point or hole) can present irregularities or de- 
fects which may induce undesirable and intrinsic changes 
in the outcoming signal. In this work we show that 
rounded borders in the open set generate a rich dynamics 
in the integrable rectangular billiard. The escape times 
statistics for the long-lived trajectories present charac- 
teristic of sticky motion when the rounded border repre- 
sents around 0.1% ^ 10% from the whole billiard size. 
In this sticky region the escape times decay exponent 
7esc increases with R/L and approaches 7esc '^2.0 for 
R/L = 1/10 {R is the radius of the rounded border and L 
is the height of the billiard) . Emission angles and escape 
times show a self-similar structure only for initial angles 
inside backgammon like stripes which are born at the in- 



tegrable case R/L = 0. Trajectories which start inside 
the backgammon stripes will collide, at least once, with 
the rounded border, generating the power law behavior. 
As R/L increases more and more, different backgammon 
stripes start to overlap and the exponent jesc ^ 2.0 is 
obtained. For R/L = 1.0 the exponential decay for the 
escape times is reached, which corresponds to the chaotic 
motion. 

From the nonlinear perspective our results are impres- 
sive, showing that a very rich dynamics and stickiness 
comes alone from tiny border effects (~ 0.1%) in in- 
tegrable billiards. Such effects, including the Arnold 
tongues, should be visible directly in experiments with 
open integrable devices, open systems with leaks (see [iq 
and therein cited references) and also in problems related 
to conduction fluctuations in semiconductors [ll|, |lj, 
particles transport in nanostructures and cold atoms in 
open optical billiards |15l |. 

Instead of rounded borders, other shapes for the open 
sets could be used. All of them should generate sticky 
motion in integrable billiards, strongly affecting the out- 
coming signal. Our results will change quantitatively 
when the aperture a increases, however, in such cases 
long-lived trajectories will disappear and the statistical 
analysis for the escape times is not adequate anymore. 
In this work we considered an integrable billiard as the 
starting point. We also performed extensive numerical si- 
mulations to study rounded border effects in the stadium 
billiard, which is already chaotic without border effects 
{R/L = 0). For all values of R/L we found that the es- 
cape times statistic has an exponential decay, expected 
for the chaotic behavior. 
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